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<N : Abstract 

^1 A unique analytical result for the Migdal-Kadanoff hierarchical lattice is obtained. The scaling 
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of the defect energy for a zero-dimensional spin glass is derived for a bond distribution that is 
continuous at the origin. The value of the "stiffness" exponent in zero dimensions, yo = — Ij 



^ \ corresponds to the value also found in one dimension. This result complements and completes 

S . 

c/3 ' earlier findings for y^ at d > 0. 
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A quantity of fundamental importance for the modeling of amorphous magnetic materials 
through spin glasses [l| is the "stiffness" exponent y (3,y|- The stiffness of a spin configuration 
describes the typical rise in magnetic energy AE due to an induced defect-interface of 
size L. In a glassy system, the potential energy function resembles a high-dimensional 
mountain landscape over its variables |j|. Any defect-induced displacement of size L in such 
a landscape may move a system equally likely up or down in energy AE. Averaging over 
many incarnations of such a system then results in a typical energy scale 

<|A^|>~L'^ (L^oo). (1) 

The importance of this exponent for small excitations in disordered spin systems has been 
discussed in many contexts [l|, |2, y, la, lO, LD, 1^- In particular, it signifies a renormalized 
coupling strength (across any hypothetical interface) between regions in space separated by 
a distance L p|: li yd > 0, regions in space are strongly coupled at low temperature and 
spin glass ordering ensues, i. e. Tg > 0. 

Ref. [si provided a description of yd as a continuous function of dimension d using a fit 
to the data obtained in Ref. (ll| for d = 2, 3, . . . , 6. That fit became credible in that it 
reproduced the exactly known result in d = 1, m = —1, to within less than 1%. Hence, 
it validated the values for yd found in Refs. [HJ, luj and produced a number of predictions 
such as that di = 5/2 may be the lower critical dimension (the dimension in which yd = 0) 
for Ising spin glasses, in accordance with an earlier calculation invoking replica symmetry 
breaking |l2l |. 

In a quest for understanding universality in spin glasses, there has been considerable 
interest recently in the behavior of yd even for d < di, where any spin glass ordering is 
unstable. Presumably, for divergent energy scales p| in Eq. ((IJ, i. e. for yd > 0, universality 
holds and low-temperature properties of the system are independent of the details of the 
bond distribution, as long as it possesses a zero mean and unit variance. In contrast, below 
the lower critical 
distributions 

leads to trivial scaling in Eq. (|H), as was found numerically iov d = 2 |l3| and is exactly 
known for d = 1. (In a linear chain of L spins, the T = energy difference \AE\ for reversed 
boundary conditions is given by the smallest bond-weight, \J\ = 1, independent of L.) Only 
bond distributions P{J) which are continuous near the origin -P(O) obtain non-trivial scaling 



dimension significant differences have been found between classes of bond 
1^, llJ, llJ|. Especially, for all d < di, a discrete bond distribution (J = ±1) 
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Figure 1: Plot of the spin glass stiffness exponent Vd as a function of dimension d. Shown are the 
values for y^, on hyper-cubic lattices from Refs. [ll|, [l^ and a cubic fit to that data. For bounded 
distributions, the exponent remains locked at y'^^^'^ = for d < di = 5/2 (upper horizontal line). 
For continuous distributions, the fit reproduces the exact result, yi = —1 to 0.8% and suggests 
yo = —2. In contrast, the result here suggests yg""* = —1 and hence, y'^""* = —1 for all d < 1 (lower 
horizontal line). 

as represented by the curve in Ref. [911, including the exact result yi = —1. (Here, the smallest 
bond- weight \J\ in the chain approaches zero with l/L.) 

In this paper we report on a (rather fortuitous) analytical result for a zero-dimensional 
spin glass that further clarifies the behavior of Ud^^ foi^ ^ continuous bond distribution. In 
the Migdal-Kadanoff hierarchical lattice (MK) |2, JS, 1^, llJ| we find for a ci = dimensional 
spin glass that yo = ^1 exactly. 

To our knowledge, aside from d = 1 (where MK is trivially exact) and the large-link 
limit |3|, this is the only exact result for MK applied to spin glasses. While not of great 
practical relevance, studying physical systems in un-physical dimensions has proved to be 
of significant theoretical relevance (iSl, |la|. A zero-dimensional spin glass in particular has 
in fact been considered previously in Ref. (l7i |. 

Our result suggests a behavior for yd as depicted in Fig. [ij While for a discrete bond 
distribution it is yf^'^ = for all d < di, for a continuous bond distribution, yd first extends 
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Figure 2: Bond-moving scheme in the Migdal-Kadanoff hierarchical lattice, here for a square lattice 
{d = 2) with / = 2, i. e. 6 = 2 in Eq. Q. Starting from the lattice with unit bonds (left), 
bonds in intervening hyper-planes are projected onto every /th plane in one direction (middle), 
then subsequent directions, to re-obtain a similar hyper-cubic lattice of bond-length I (right). The 
renormalized bonds in this case consist of 6 = 2 branches, each of a series of / = 2 bonds. 

smoothly to (non-trivial) negative values through di tovv^ards d = 1, beyond vi^hich it appears 
to get fixed at yf^^ = —1 for all c? < 1. As long as y^ > for d > di the exponent is believed 
to be universal, independent of the bond distribution. 

The Migdal-Kadanoff (MK) hierarchical lattice (l^,l20| provides a real-space renormaliza- 
tion scheme that approximates especially low-dimensional spin glasses well and is trivially 
exact m d = 1. These lattices have a simple recursive, yet geometric, structure and are 
well-studied |2, y, la, Isjj . Starting from generation I with a single bond, at each subsequent 
generation / -|- 1, all bonds from / are replace with a new sub-graph. This structure of the 
sub-graph arises from the bond-moving scheme, as shown in Fig.[2l in d dimensions |l9l.l2ril|: 
In a hyper-cubic lattice of unit bond length, at first all / — 1 intervening hyper-planes of 
bonds, transverse to a chosen direction, are projected into every Ith hyper-plane, followed by 
the same step for / — 1 hyper-planes being projected onto the /th plane in the next direction, 
and so on. In the end, one obtains a renormalized hyper-cubic lattice (of bond length /) in 
generation / -|- 1 with a reformulated (/ -|- l)-bond consisting of a sub-graph of 
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parallel branches of a series of / /-bonds each. We can rewrite Eq. 

ln(6) 



as 



d= 1 + 



HI) 



(2) 



(3) 



anticipating analytic continuation in / and b to obtain results in arbitrary dimensions d. 

Instead of solving the problem on the hyper-cubic lattice, we merely need to consider the 
recursive scheme of obtaining the bond-distribution in generation / + 1 from sub-graphs of 
bonds from generation /. Numerically, this is done efficiently at any temperature T to yield 
a stationary bond distribution for J — > oo, i. e. the thermodynamic limit L = /^ — > cxd. 

Here, we are only concerned with T = 0, which simplifies the calculation drastically to 
the point that analytical results can be obtained. For instance, a series of / bonds can always 
be replace by the bond of smallest absolute weight. Thus, if these bonds are drawn from a 
distribution Pi (J), then the distribution Qj (K) of the effective bond K replacing the series 
can be obtained from 

Q^PiK) =V{K = sign (Ji x J2 x . . . x J,) ™« (l^il , • • • , \Ji\)] , 

/OO /"OO 

dJiPi (Ji) / dJi^^Pj ( J,_i) e i\Ji\ - I J,_i|) ... (4) 
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dJiPi{Ji)Q{\J2\-\Ji\)5{\Ji\-\K\), 



where Q{x) refers to the unit step-function and 5{x) = <d'{x) is the Dirac delta- function. In 
writing Eq. (^, we have already exploited the symmetry of the integrand under relabeling 
Ji, . . . , J/, and we have dropped the obvious norm of Q^^\K). 

Similarly, having h bonds drawn from a distribution Q^^\K) in parallel leads to the 
distribution P/+i(J) of the reformulated bond at generation / + 1, 

/OO '' / b ' 
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where we have used the integral representation 6{x) = f^^due'^'^^ /{2tt). 

For general Pj, Eqs. (JUIHI) both are quite complicated. This is particularly true, since we 
would be interested in finding the limiting shape of Pi{J) after infinitely many iterations, 
/ — > 00. In that limit, we expect the effective coupling between domains a distance L = f 
apart to scale as p| 

Je//(L) ~ l^JefAL/l). (6) 

Each iteration the characteristic "width" Jeff{L) = [(J^")^]^/^™ ^ {\AE\)^, of Pi{J) in- 



creases (or decreases) by a factor of l^, so 

Pi+iiJ) 
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Clearly, if the width grew smaller, /^ < 1, the behavior of Pi{J) near J = would become 
increasingly relevant, explaining the non-universal behavior there ^]. 

We can now report on a solution for Eqs. (jHEI) for all values of / for a specific set of 
initial distributions and choice of / and 6, appropriately analytically continued. With those 
choices, only one iteration of Eqs. (jMSI) is necessary, since the distribution remains shape 
invariant, i. e. Eq. becomes an equality. This set of distributions is continuous and finite 
near the origin and thus should represent the universality class containing Gaussian bonds, 
for instance. Starting more generally with an initial distribution 
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where Jo > sets the energy scale for this distribution, it is easy to show that 
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by recursion of Eq. (jij). Note, that Eq. readily continues to any real value of / > 0. 

The evaluation of Eq. ^ is somewhat more complex. To facilitate the subsequent anal- 
ysis, it is best to consider the moment generating function for Pi (J), 
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where the integral over J merely represents a 5— function. Using Eq. (jH)), we find 

POD "I b 

h[a) = I dKQf{K)e-''''' 



/ (g + 1) / dx{l~x)'^^ cos (a Jqx) 
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(11) 



Correspondingly, we find for the generating function of the initial bond distribution in 

Eq. JED, 

• 1 



,(«) = (e-'"-^)^ ,= {q + l) dx {1- xy cos (a Jqx) . 
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Despite the obvious similarities between Eqs. ()11II12|) . finding a set of parameters that 
make 0o and 0i similar are hard to find because of the exponent b in Eq. (fTT| . Even for 
6=1 and I ^ 1, i. e. d = 1 according to Eq. Q, we need to iterate for / ^ oo to find the 
width to decay to zero with 1/L, i. e. yi = —1. It is merely a lucky circumstance that there 
are two solution with the required properties, each independently yielding the same result. 

First, for / = 2 we have to continue the "branching number" to 6 = 1/2 to obtain d = 
in Eq. ^. Then, for the rectangular function for Po{J), i. e. q = 0, we find 



sin (g Jq) 

(/)o(a) = , (pi [a] 



4sin2(^) 



(aJo)^ 



(13) 



both of which are invariant in an sufficiently large open interval around a = 0, required to 
generate any moment. We identifying in Eq. (P) Jeff{L = /°) = Jq and J'eff{L = l^) = 
Jo/2 = PJeffiL = f), hence, yo = -I- 

Second, for 6 = 2 we have to continue the "series number" to / = 1/2 to obtain ci = in 
Eq. Q. Then, for the triangular function for Po{J), i. e. q = 1, we find 

^ ^ ^ 4sin2(s^) sin^aJo) 

[aJo) [aJo) 

In this case, we identify in Eq. (P) again J'eff{L = P) = Jq, but Jeff{L = l^) = 2Jo, which 
now results again in yo = —1 because length-scales are actually shrinking, / = 1/2. 

It is clear that nothing will change on this result under further iteration, / — i> / + 1. 
Assuming that there is a unique solution for Eq. ^, one would expect that any choice for 
Po{J) that is continuous at J = should converge to a rectangular (triangular) function for 
/ = 1/6 = 2 (= I). We have not been able to extract any further result of this nature. 

Hence, we are lead to believe that yo = —1. Of course, such a result has to be considered 
with care. It is a considerable weakness of MK that its results for a given dimension d 
according to Eq. ^ are not generally unique [3|. Here, at least, we found two different 
combinations of / and 6, both giving d = and yielding an identical stiffness exponent. 
This may indicate a unique result for any combination / = 1/6. Furthermore, it is not clear 
that the MK result for c? = should necessarily correspond quantitatively to an Edwards- 
Anderson spin glass on a zero-dimensional lattice. This is only known to be true iov d = 1 
(i. e. 6=1) and certainly wrong for any d > 1, albeit with slowly increasing error. 
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